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We examine the general features of the non-commutativity of the magnetization operator and 
Hamiltonian for the small quantum spin clusters. The source of this non-commutativity can be a 
difference in the Lande g-factors for different spins in the cluster, XY-anisotropy in the exchange 
interaction and the presence of the Dzyaloshinskii-Moriya term in the direction different from the 
direction of the magnetic field. As a result, zero-temperature magnetization curves for small spin 
clusters mimic those for the macroscopic systems with the band(s) of magnetic excitations, i.e. for 
the given eigenstate of the spin cluster the corresponding magnetic moment can be an explicit func¬ 
tion of the external magnetic field yielding the non-constant (non-plateau) form of the magnetization 
curve within the given eigenstate. In addition, the XY-anisotropy makes the saturated magnetiza¬ 
tion (the eigenstate when all spins in cluster are aligned along the magnetic field) inaccessible for 
finite magnetic field magnitude (asymptotical saturation). We demonstrate all these features on 
three examples: spin-1/2 dimer, mixed spin-(l/2,l) dimer, spin-1/2 ring trimer. We consider also 
the simplest Ising-Heisenberg chain, the Ising-XYZ diamond chain with four different g-factors. In 
the chain model the magnetization curve has a more complicated and non trivial structure which 
that for clusters. 

PACS numbers: 75.10.Pq, 75.50.Xx 

Keywords: low-dimensional quantum magnetism, magnetization plateaus, molecular magnets, magnetic 
anisotropy 


I. INTRODUCTION 

Magnetization curves of low-dimensional quantum an- 
tiferromagnets are topical issue of current research inter¬ 
est, because they often involve intriguing features such as 
magnetization plateaus, jumps, ramps and/or kinks. The 
spin-1/2 quantum Heisenberg chain, the spin-1/2 quan¬ 
tum Ising chain in a transverse field, and the spin-1/2 
quantum XX chain in a transverse field are a few paradig¬ 
matic examples of exactly solved quantum spin chains for 
which zero-temperature magnetization varies smoothly 
with rising magnetic field until the saturation magne¬ 
tization is reachedji^— Contrary to this, the integer- 
value quantum Heisenberg chains (and also many other 
low-dimensional quantum antiferromagnets) contain in 
a zero-temperature magnetization process remarkable 
magnetization plateau(s) at rational value(s) of the sat¬ 
uration magnetizationi^i^ The intermediate plateaus of 
Heisenberg spin chains reflect quantum states of mat¬ 
ter with exotic topological order such as the Haldane 
phasej^ii whereas their presence is restricted by quan¬ 
tization condition known as Oshikawa-Yamanaka-Affleck 
rule.— 

On the other hand, it could be generally expected that 
the antiferromagnetic Heisenberg spin clusters should al¬ 
ways exhibit leastwise one intermediate plateau before 
the magnetization jumps to its saturation valuej^^— 
This naive expectation follows from the energy spectrum 
of the quantum Heisenberg spin clusters, which is com¬ 


posed of a few discrete energy levels that cannot natu¬ 
rally form a continuous energy band needed for a smooth 
variation of the magnetization at zero temperature. At 
first sight, this argumentation is consistent with the ex¬ 
istence of at least one plateau and magnetization jump, 
which bears a close relation to level crossing caused by 
the external magnetic held. From this perspective, the 
quite natural question arises whether or not intermediate 
magnetization plateau(s) can be partially or completely 
lifted from zero-temperature magnetization curves of the 
Heisenberg spin clusters. 

Another spin systems which should be noted in the 
context of the small quantum spin clusters are the Ising- 
Heisenberg chains. They are the one-dimensional spin 
systems where the small quantum spin clusters are as¬ 
sembled to the chain by alternating with the Ising spins 
in such a way, that the Hamiltonian for the whole sys¬ 
tem is a sum of mutually commuting block Hamiltonians. 
These systems have much in common with the “classical” 
chains of the Ising spins, as they can be solved by the 
same technique and the eigenstates are just the direct 
product of the eigenstates of the single block, though, 
for the more complicated structure the doubling of unit 
cell is possible. The magnetization curves, thus, for the 
Ising-Heisenberg spin systems share almost all features 
with the magnetization curves of the small spin clus¬ 
ters, but can contain much more intermediate magneti¬ 
zation plateaus. Various variants of the Ising-Heisenberg 
chains have been examined: diamond-chaiiiji^^— saw- 
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tooth chain )^^i^° orthogonal-dimer chain^^— tetrahe¬ 
dral chain)^^— and some special examples relevant to 
real magnetic materials.—^— 

In the present work, we will rigorously examine a mag¬ 
netization process of a few quantum Heisenberg spin clus¬ 
ters and Ising-Heisenberg diamond chain, which will not 
display strict magnetization plateaus on assumption that 
some constituent spins have different Lande g-factors 
and may be a XY-anisotropy of the exchange interac¬ 
tion. Also, the Dzyaloshinskii-Moriya (DM) term in a 
direction different from that of the magnetic field can 
lead to the same effect. All those features of the spin 
Hamiltonian make the magnetization non-conserved, i.e. 
non-commutating with the Hamiltonian. This specific 
requirement naturally leads to a nonlinear dependence of 
the energy levels on a magnetic field, which consequently 
causes a smooth change of the magnetization with the 
magnetic field within one and the same eigenstate. Al¬ 
though the smooth change of magnetization due to a 
difference in Lande g-factors or/and XY-anisotropy and 
non-collinear DM-term may be quite reminiscent of that 
of quantum spin chains with continuous energy bands, it 
is of course of completely different mechanism with much 
simpler origin. 

The single-chain magnet, [{(CuL) 2 Dy}{Mo(CN) 8 }] • 
2 CH 3 CN • H 20 ^“— is a remarkable example of both 
Ising-Heisenberg one-dimensional spin system and a spin 
model with different Lande g-factors, leading to non¬ 
plateau form of the region of the magnetization curve 
corresponding to the same eigenstate. However, as the 
exact analysis shows^ the effect is just hardly visible 
in the magnetization curve plot in virtue of the very 
small difference in Lande g-factors of the magnetic ions, 
though, the exact expression for the magnetization has 
explicit dependence on the magnetic field. Almost the 
same effect but even quantitatively less pronounced have 
been observed in the approximate model of the one¬ 
dimensional magnet, the F-F-AF-AF spin chain com¬ 
pound Cu( 3 -Chloropyridine) 2 (N 3 ) 2 i^ 

The organization of this paper is as follows. In the 
next Section, we will clarify a few general statements 
closely related to an absence of actual plateaus in zero- 
temperature magnetization curves of quantum spin clus¬ 
ters and chains. These arguments of general validity will 
be subsequently illustrated on a few specific examples of 
the spin-1/2 quantum Heisenberg dimer, the mixed spin- 
(1/2,1) Heisenberg dimer, the spin-1/2 Heisenberg trimer 
and the spin-1/2 Ising-Heisenberg diamond chain in the 
following four Sections. The summary of the most impor¬ 
tant findings along with the implications for experimental 
systems will be presented in the concluding part. 


II. GENERAL STATEMENTS 

Let us first start with a few very general statements 
elucidating the issue of the non-constant magnetization 
within one physical state or the explicit magnetic held 


dependence of the magnetization corresponding to a cer¬ 
tain eigenstate of the small spin clusters. Obviously, the 
aforementioned phenomenon arises when the magnetiza¬ 
tion is not a good quantum number 

[n, M^] ^ 0 , ( 1 ) 

Here, 'H. stands for the Hamiltonian of a spin cluster. One 
can distinguish two cases, when z-projection of the total 
spin S'/qj does not commute with the Hamiltonian and 
the magnetization operator is proportional to it 

( 2 ) 

or when the z-projection of the total spin is a good 
quantum number, but the magnetization operator is not 
proportional to it and does not commute with the Hamil¬ 
tonian 

[■^5 'S'tot] = 0; 7^ 9^J'BStot■ (3) 

Of course, another possibility is to have the magne¬ 
tization which is non proportional to and the z- 
projection of the total spin non-conserved. The spin 
Hamiltonians, which do not commute with usually 
contain XY-anisotropy or/and DM vector with a non¬ 
zero X or Y part. The magnetization is non-proportional 
to the total spin when the spins possess different 
Lande g-factors. 

III. SPIN-1/2 HEISENBERG DIMER 

In this section we consider the spin-1/2 Heisenberg 
dimer as the simplest system of two interacting quantum 
spins described by the most general Hamiltonian 

-Harm = J{{l + i)SlS^ + [l--f)SlSl + ^SlS^^} 
+D • (Si X S 2 ) - B • (51 Si + 32 ^ 2 ). (4) 

Here, S“ 2 i (a = x,y,z) are the spatial components of 
the spin-1/2 operators for two spins in the dimer. We 
assume the fully anisotropic XYZ Heisenberg coupling 
with two anisotropy constants 7 , A and two different but 
isotropic Lande g-factors. The spatial direction of the 
magnetic field B and the DM-vector D are arbitrary so 
far. Without loss of generality, one may however choose 
a direction of the magnetic field along the z-axis and the 
DM vector to lie in xz-plane 

nd^m = J{{i + ^)s^s^ + {i--f)sfsi + As^,si} 
+D, {SfS^ - S/S|) + D, {S^S^ - SfS^) 
-B{g^S! + g2SI). (5) 

Let us calculate the commutators of the Hamiltonian ([S]) 
with the z-projections of the operators corresponding to 
the total spin and magnetization 

SL = S! + S^, ( 6 ) 

= giSl+g2SI, 
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[Hd.m, = - 2 z7 (S^S^ + SfS^) + iD^ (S^S^ - S^S^), 
in^^rn, M^] = - ig-D, ( 5 J= 52 " + SlSl) + 

+ ig.J{SlSl-S\S^^) + 

-ng+ {S^Sl + SfS^) + 

+ iD^{g,S^,S^-g 2 S^,S^), ( 7 ) 

where g± = gi± g2. As one can see, the XY-anisotropy, 

7, and the DM-vector x-projection, make the S^g^ 
and Ai^ non conserved, but even if we set them to zero, 
the magnetization may still be a non-conserved quantity 
because of the difference in Lande g-factors. Thus, the 
spin- 1/2 Heisenberg dimer may exhibit the non-constant 
magnetization within one ground state if at least one of 
the parameters, 52 — ffi, 7 or Dx{Dy) is non-zero. Let us 
put Dx = 0 as it makes the analytic calculations quite 
cumbersome (the eigenvalue problem leads to a quartic 
equation), and start with the exact diagonalization of 
the Hamiltonian for the anisotropic spin- 1/2 Heisenberg 
dimer with different Lande g-factors. The eigenvalues are 

£ 1,2 = -^ ± yB^gl + J^ + Dl 

e,^, = y±yB^gl + J^^A (8) 

The corresponding eigenvectors are 


1^1,2) = 




(I ti) + A±| 4,t)) : 


A± =p±e^'^, (j) = arctan 


£1 

J 


P± =- 


Bg_±JB‘^gl + J^+Dl 


1^3,4) = 


1 + Bl 


+ Dl 

(ltt) + i?±IU)): 


B± = 


Bg+ ± ^B^gl + 


J7 


( 9 ) 


Under the conditions gi = g2 (and = 0 ) the first 
two eigenstates become conventional singlet and = 0 
component of the triplet respectively. However, there is 
no continuous transition to the = 1 and = — 1 

components of the triplet in at 7 —>■ 0 . In order 

to obtain | ft) a-nd | J,),) one has to put 7 = 0 in the 
Hamiltonian before diagonalization. Let us calculate the 
magnetization eigenvalues for all those eigenstates: 

(^i,2i (51^/+ 52 ^ 1 ) i'i'i,2)= 

2 1 +Pi, 

-I- 71 n^ 

“ ^2^B2g2 _,_j2+£,2 ■ y ^1 

At gi = g2 this expression becomes 0 . However, at 
gi 7^ 52 we have explicit dependence of the eigenvalue, 



B 

FIG. 1 : The energy spectrum of the isolated S'=l /2 dimer 
with gi = 2 , 32 = 3 , J = 1 , Dj = 1 , 7 = 2 and A = 3 
displaying level crossing. Two bottom curves corresponds to 
£2 and £4. The non-linearity in B on the energy levels is the 
main reason for the non-plateau magnetization. 


- y= 0 y=l y=2 y=4 



FIG. 2 : The zero temperature magnetization curves for the 
5 = 1/2 dimer with 31 = 2 , 32 = 6, J = 1 , Dz = 1 , A = 2 and 
7 = 0 (red solid); 7 = 2 (black dotted); 7 = 4 (blue dashed) 
and 7 = 6 (brown dot-dashed). Maat = ^(31+32) = 4 . 


corresponding to the certain eigenstate on the magnetic 
field. This leads to a non constant magnetization for 
the given eigenstate. Thus, we have here = 0 and 
yf 0. 

For the other two eigenstates we have 


('^'3,4! {giSi + g 2 Si) I'I'st) — 2^'+ 1 _|_ 52 


= T 


Bgi 




( 11 ) 


At 7 = 0 the expression transform to 


('1'3,4| {giSi + 32<S'2) — Y25+) 


( 12 ) 
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which corresponds to | ft) and | J'J-) eigenstates. The 
Eq. (fTTl) has another important feature. The transverse 
quantum fluctuations enhanced by the XY-anisotropy 7 
reduce the magnetization in z-direction in such a way 
that it never reaches its saturated values ±5(51 +32) at 
any non-zero 7 and finite magnetic field B. It is also 
important that even at the equal g-factors, 02 = gi = 9 
the magnetization expectation values for the eigenstates 
|4>37) exhibit explicit magnetic field dependence and do 
not reach their saturated values at non-zero 7. Another 
case of interest is the 92 = —91 = 9, when the magne¬ 
tization expectation value for the I'I'i 2) is non-zero and 
exhibits explicit dependence on the magnetic filed and 
the eigenstates In's,4) demonstrate zero magnetization. 
Moreover, the corresponding expectation values become 
singular at 7 = 0, because, as was mentioned above, there 
is no continuous limit 7 —>■ 0 in terms of eigenvalues and 
eigenvectors. The magnetic susceptibility for the afore¬ 
mentioned states can be obtained in a straightforward 
way by taking a derivative of the Eqs. m and m 
with respect to B: 


^(Vl/3.4|(5l5r+52^|)|'f'3,4) 




J" 7 ^g+ 

'2^(324+72^2)3- 


( 13 ) 


The deviation form the horizontal line for the zero- 
temperature magnetization curve of the dimer under con¬ 
sideration, thus, is governed by three factors. Difference 
of Lande g-factors, 92 — gi is responsible for the non¬ 
plateau behavior for the initial part of the magnetization 
curve, B < Be, the larger is the absolute value of the 
difference the more pronounced the deviation is. At the 
same time, the overall Lande g-factor, gi + 92 and the 
XY-anisotropy 7 make another part of the magnetiza¬ 
tion curve, which in the limit 7 = 0 correspond to the 
saturation, | tt)j non flat. The critical field B^ is found 
from the level crossing. As for the zero temperature only 
|4'2) and |'I'4) are realized we can found the correspond¬ 
ing value of the magnetic field from the equation, £2 = £4, 
which leads to 


J 


Be 


^2gig2T- - 1 - JA^{gf+g^)JA + \J^g\gl {{gf+g^)r- +2gig2T+) + {gf - g^)"^ 

‘^V^9i92 


( 14 ) 


with 


r± = (D 2 + j 2 (l± 72 )). ( 15 ) 

The typical picture of the level crossing curve one can 
see in Fig. [H However, the B = 0 ground state is also af¬ 
fected by the value of the XY-anisotropy 7. The ground 
state becomes |4'4) for a sufficiently large 7 above cer¬ 
tain critical value je- The critical value is given by the 
equation 


7e = A + y'l + {D,/J)^. ( 16 ) 

The non-linear behavior with respect to the magnetic 
field is the main reason for the non-plateau magnetiza¬ 
tion. As the DM-term in z-direction does not bring any 
qualitatively new physics, we can hereafter put = 0. 
The expression for the critical field (fl 3 does not lead to 
a proper 7 = 0 limit. The case of isotropic Heisenberg 
interaction must be considered separately. In the case 
of isotropic Heisenberg interaction the difference is only 


in the saturated states presented here, and which trans¬ 
forms to |4>37) at non-zero 7, while the |4>i_2) eigenstates 
remain the same. The value of critical field in this case 
is 


5 +A -f 
Be = J - 


glA'^ + 49192 
4 gi 52 


( 17 ) 


Thus, the jump to the saturated magnetization takes 
place for the 7 = 0 at this value of the magnetic field. 
The magnitude of the jump depends on the difference of 
the Lande g-factors and is given by 


AM = ^ll- 


4 gig 2 \/1-1- 


9 ^ 

- 



3 + 



9^ [s+A+^g^ ^2+49192] 

issfgl 


}■ ( 18 ) 


The corresponding plots of the zero-temperature mag¬ 
netization one can find in the Figs. [ 2 ] and [H In the 
Fig. [ 2 ] the evolution of the T = 0 ground state for dif¬ 
ferent values of 7 are presented for J = 1 , = 1 , 
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A = 2 , 5i = 2 and g2 = 6. The critical value of 7 at 
which the B = 0 ground state of the spin- 1/2 spin dimer 
changes from |'I'2) to |d>4) for these values of J,Dz and 
A is 7c = 2 -I- y /2 ~ 3 . 41 . Therefore, for 7 = 0 and 
7 = 2 one can see magnetization curves with two eigen¬ 
states separated by the jump. The non-plateau behavior 
of the magnetization for |4'2) at B < B^ is well visi¬ 
ble. Also, the non-plateau character of the magnetization 
curve, corresponding to |'I'4) is obvious for 7 = 2, while 
for 7 = 0 we see ideal plateau M = ^[gi + 52) = 4 . 
This non-plateau behavior and inaccessibility of the sat¬ 
uration are more pronounced for 7 = 4 and 7 = 6 when 
the system for all values of the magnetic field {B > 0 ) 
is in the |'I'4) eigenstate and its T = 0 magnetization 
curve demonstrates the form very similar to the one of 
a system with a band of magnetic excitations or/and to 
the high-temperature magnetization curve given by the 
Brillouin function. The effect of the difference between 
the Lande g-factor is summarized in Fig. O To demon¬ 
strate the evolution of the ground state |'I'2) under the 
change of the difference of Lande g-factors, we have cho¬ 
sen 7 = 0 and J = 1 , = 0 , A = 1 and plotted the 

normalized magnetization M/Mgat-, as saturation magne¬ 
tization, Msat = ^(,91+92) is different for each curve. For 
the gi = g2 there are just two ideal plateaus at M = 0 
(singlet state) and M = 1 . The magnetization jumps 
from 0 to 1 at Be = (here 9 i = 92 = g)- How¬ 

ever, the essential changes appear when the difference 
between g-factors is growing. For 52 — gi non equal to 
zero the part of the magnetization curve corresponding 
to the -I- = 0 deviates from the horizontal line 

and become almost linear (for small g_) and then more 
and more rapidly growing with the shift of the transition 
point between |4'2) and |4'4) in the lower B region. 

Of course, the effects of the DM-terms in molecular 
magnets and other low-dimensional many-body spin sys¬ 
tems have been intensively studied in various contexts 
during last decade 42 i“— In Refidl the isolated spin-dimer 
with DM-terms has been considered with general mutual 
orientation of the DM-vector and magnetic field. 


B = J Si • S2 + D • (Si X S2) - gB{S^ + B|), ( 19 ) 

where D = { 0 , D sm 9 , D cos 9 ). Though, as was shown 
above, in this case the eigenvalue problem leads to the 
solution of quartic equation, the authors found an ap¬ 
proximate ground state in the limit D/J -C 1 and below 
the critical field B = J/g. They explicitly found out the 
magnetization of the ground state which turned out to 
be linear in B 


It is also straightforward to construct the thermody¬ 
namics of the simple dimer. The partition function is 
calculated directly from the spectrum: 


Zdim =2 ch 

+ e ^ ^ ch 


^B^gl + J^ + Dl 


^^B^gl + JW 


( 21 ) 


The magnetization is found in a standard way, as Mdim = 
1 yielding 


^dim — 



r 4-^- .u 

/3^S2g2_^J2^2 

Zdtm 


2 


^B^gl+J^+Dl 


sh 


0y/B^^T+J^TD% 

2 


■ ( 22 ) 


The plots of the finite-temperature magnetization for 
the S=l /2 spin dimer are presented in Fig. | 4 l It is 
worth mentioning that thermal fluctuations eliminate 
from magnetization curves all typical structures (such as 
plateaus or quasi-plateaus) quite similarly as the large 
XY-anisotropy does for zero-temperature magnetization 
curves. 


IV. MIXED-SPIN HEISENBERG DIMER 

Another interesting example of a simple quantum 
spin system, which may possibly show a striking depen¬ 
dence of the total magnetization on a magnetic field, is 




B/J 


M = ^ (D X B) X D. ( 20 ) 

This is the approximate form of the non-linear behavior 
of the magnetization we have obtained exactly above in 
the case of D = ( 0 , 0 , D^). Despite of all these results, the 
issue of the non-conserving magnetization and its conse¬ 
quences have not been systematically investigated so far. 


FIG. 3 : The normalized zero-temperature magnetization 
curves, M/Maat of the S = 1/2 dimer with two different 
g-factors in case of isotropic exchange interaction, 7 = 0 
and A = 1 . Here, for the sake of simplicity, we put 
J = l,Dz = 0 ,gi = 2 and present the curves for the dif¬ 
ferent values of 52- From the bottom to top g2=2 (red); 4 
(green); 6 (blue); 8 (magenta); 10 (brown) and 20 (black). 
Msat = 5(51 +52). 
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FIG. 4: The normalized finite-temperature magnetization 
curves, M/Msat for the S' = 1/2 spin dimer with two different 
p-factors in case of isotropic exchange interaction, 7 = 0 and 
A = 1 at different temperatures for <71 = 2 , 32 = 6 , J = l, 
= 7 = 0. T/J = 1.7 (red, solid); T/J = 0.1 (black, 
dotted); TjJ = 0.03 (blue, dashed) and T/J = 0.001 (brown, 
dot-dashed). Msat = \{gi + 92 ) =4. 


the mixed spin-(l/ 2 ,l) Heisenberg dimer defined by the 
Hamiltonian 


+ -B{g,Sl+g2^il). (23) 


Here, Sf and /xf {a = x, y, z) represent spatial compo¬ 
nents of the spin-1/2 and spin-1 operators, respectively, 
the exchange constant J denotes the XXZ Heisenberg 
coupling between the spin-1/2 and spin-1 magnetic ions, 
A is an exchange anisotropy in this interaction, D is an 
uniaxial single-ion anisotropy acting on a spin-1 magnetic 
ion, gi and g2 are Lande 5-factors of the spin- 1/2 and 
spin-1 magnetic ions in an external magnetic field B. As 
the effect of the XY-anisotropy was described in details 
in the previous Section, here to put it simple, we assume 
7 = 0 . A straightforward diagonalization of the Hamilto¬ 
nian (1^ gives a full spectrum of eigenstates, which can 
be characterized by the following set of eigenvalues 


and the corresponding eigenvectors 

|«'i,2) = iTiTl), 

l'f'3.4) = c±|-i,l)Tcf|i,0), 

I'I's.e) = c^l~1) T cjl ~ 0), (25) 


where the respective probability amplitudes are given by 


Cl = 


ct = 


1± 


JA-2D-2g-B 


^(JA- 2 n- 25 _B) 2 -|- 8 J 2 

/I ^ ,JA-2D+2g^B 

2 ^(JA-2n+2g„B)2-|-8J2 


( 26 ) 


It should be mentioned that the mixed spin-( 1 / 2 , 1 ) 
Heisenberg dimer exhibits a strict intermediate plateau 
at one-third of the saturation magnetization regardless 
of uniaxial single-ion anisotropy on assumption that the 
Lande g-factors of both constituent magnetic ions are 
equal 51 = 52- If the Lande g-factors are different 
gi ¥^92, then, the mixed spin-( 1 / 2 , 1 ) Heisenberg dimer 
displays more intriguing zero-temperature magnetization 
curves basically affected by a relative strength of the uni¬ 
axial single-ion anisotropy. To illustrate the case, the to¬ 
tal normalized magnetization of the mixed spin-(1/2,1) 
dimer is plotted in Fig. [ 5 ] against the magnetic field for 
the isotropic Heisenberg coupling A = 1 , several values 
of the uniaxial single-ion anisotropy and two particular 
sets of Lande g-factors. A smooth variation of the total 
magnetization observed below a saturation field relates 
to a gradual change of probability amplitudes of two en¬ 
tangled microstates | — 1/2,1) and |l/2,0) within the 
eigenstate l^'a). In a low-field region with continuously 
varying magnetization, mean values of two constituent 
spins and the total magnetization can be therefore cal¬ 
culated with the help of the corresponding lowest-energy 
eigenvector IT's) given by Eq. ( 1251 ) 




1 


JA-2D-2g-.B 


2 y/{JA-2D-2g-B)^+8B ’ 




1 


1 + 


JA-2D-2g^B 


y/(JA-2D-2g-BY+8J'2 


M = + g 2 g¥)\^¥ = 

^ _ 9 - JA-2D-2g-B 

2 2 y/(JA-2D-2g-B)'2+8J2 ' 


1 B 

£1.2 = - JA-1-£) =F Y (51 + 252 ), 

£ 3,4 = — — (JA — 2Z1-|- 252 I?) 

-2D- 2g-Bf + 8.P, 

£ 5,6 = —-{JA — 2D — 2g2B) 

T77/(JA -2D + 2g_Bf -b 8J2 (24) 


At 5_ = 0 we have here M = | (51 = 52 = ff)- It can 
be seen from Fig. [ 5 ] that the respective field variations 
of the total magnetization (normalized by the saturation 
magnetization Msat = 5(51 + 252)) depend basically on 
whether the Lande g-factor of the spin -1 magnetic ion 
is greater or smaller than the g-factor of the spin- 1/2 
magnetic ion. The total magnetization is gradually sup¬ 
pressed by an increase in the single-ion anisotropy in the 
former case gi < 52 (see Fig. [ 5 ^), while the total mag¬ 
netization is conversely enhanced by an increase in the 
single-ion anisotropy in the latter case 51 > 52 (see Fig. 
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FIG. 5: Zero-temperature normalized magnetization curves 
of the mixed spin-(l/2,l) dimer by assuming the isotropic 
Heisenberg coupling A = 1, several values of the uniaxial 
single-ion anisotropy and two different sets of the Lande g- 
factors: (a) gx = 2, 32 = 4; (b) gi = 4, 32 = 2. Msat = 
\{gi + 252 ). 


[S}d). In general, the total magnetization displays a con¬ 
siderable dependence on a magnetic field for small enough 
single-ion anisotropies D/J k. 0, while one recovers a 
quasi-plateau dependence with only a subtle variation of 
the total magnetization in two limiting cases D/J —^ ±00 
at which the following asymptotic values are reached 


lim M/Msat 
D / J—^co 

lim M/Msat 
D / J—¥ — (yo 


9i 

2ff2 + gi ’ 
2g2 - gi 
2g2 + gi ’ 


(27) 


quantum entanglement of two microstates |l/2, — 1) and 
I —1/2, 0) may eventually become the lowest-energy eigen¬ 
state with a positive value of the total magnetization in 
spite of negative value of the total spin = —1/2. It 
is quite evident that a strong enough easy-axis single¬ 
ion anisotropy suppresses the occurrence probability of 
the microstate | — 1/2,0), whereas the other microstate 
|l/2,—1) may lead to a positive magnetization due to 
much greater the Lande g-factor of the spin-1/2 mag¬ 
netic ion gi ^ g 2 than that of the spin-1 magnetic ion. 
Mean values of two constitutent spins and the total mag¬ 
netization follow from the corresponding lowest-energy 
eigenvector I^Fs) given by Eq. (1251) 


(«' 5 |M 2 I^ 5 ) 


1 

2 


1 

2 


JA-2D+2g^B 
JA- 2 D-|- 2 g_B) 2 + 8 J 2 ’ 

_ JA-2D+2g^B 

Y^(JA-2r>+2g_S)2-|-8J2 


M = (giS'/ -b g 2 M 2 ) l^'s) = 

^ 9- JA-2D-2g^B 

2 2 i/(JA-2n-2g_B)2+8J2 ' 


(28) 

(29) 


In the analogy with the Idta) case, here we get M = —^ 
at 9i = 92 = 9- The zero-temperature magnetization 
curves displayed in Fig. |6] afford a convincing proof that 
the total magnetization may vary continuously in a low- 
field region, then it may show an abrupt jump to an 
intermediate-field region with another continuously vary¬ 
ing magnetization terminating just at the saturation field 
(see the magnetization curves for D/ J = —0.4 and —0.5). 
In accordance with the previous argumentation, the to¬ 
tal magnetization follows the formula (l^5l) in the low-field 
region attributable to the eigenstate Idts), while it varies 
according to Eq. (I29D in the intermediate-field region 
attributable to the eigenstate Idta). The magnetization 
part corresponding to the eigenstate Idta) gradually di¬ 
minishes as the easy-axis single-ion anisotropy strength¬ 
ens (i.e. it becomes more negative) and hence, the total 
magnetization shows below a saturation field only a sin¬ 
gle region with continuously varying magnetization due 
to the striking lowest-energy eigenstate |4>3) with a neg¬ 
ative total spin but a positive total magnetization (see 
the magnetization curves for D/J = —0.6 and —1.0). It 
is straightforward to calculate the susceptibility for the 
eigenstates l^'a), and I'I's): 


[giS/ + 929^2) l ' f ' 3 , 5 ) 

^ ^2^2 {J^- 2 D- 2 g.Bf + 4 


(30) 


(JA -2D- 2g_By + SP 


1 3/2 ’ 


However, the most surprising zero-temperature depen¬ 
dence of the total magnetization can be found when the 
g-factor of the spin-1/2 magnetic ion is much greater than 
the g-factor of the spin-1 magnetic ion (gi ^ g 2 ) and the 
uniaxial single-ion anisotropy is of easy-axis type D < 0. 
It turns out that the eigenvector jit's) characterized by a 


which become zero only when gi = g 2 ov D/J —^ 00 
or A —>■ 00 . It has been demonstrated that a smooth 
variation of the total magnetization at zero temperature 
within one and the same eigenstate requires a differ¬ 
ence between the Lande g-factors. From this perspec¬ 
tive, our theoretical predictions could be more easily ex- 
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The Hamiltonian can be diagonalized in a straightfor¬ 
ward way. The eigenvalues are: 


£ 1.2 = ^ T 2 (51 + 252) B, 

£3.4 = —^ (2 -I- A) =F ’ 

£5.6 = ^ (1 - A) ± Q+ -I- 2^2-S, 
£7.8 = ^ (1 - A) ± g_ - ^g2B, 

where 

g± = i\/2j2 + (j± 5 _i?f. 

The eigenvectors are 


( 32 ) 


( 33 ) 


FIG. 6: Zero-temperature normlized magnetization curves 
of the mixed spin-(l/2,l) dimer by assuming the isotropic 
Heisenberg coupling A = 1, several values of the uniaxial 
single-ion anisotropy and the Lande g-factors gi = 6 and 
g2 = 2. Msat = |((?i + 2 ^ 2 ). 


perimentally tested for the mixed spin-( 1 / 2 , 1 ) Heisen¬ 
berg dimer, which represents plausible model for het- 
erobimetallic dinuclear complexes naturally having two 
unequal Lande g-factors due to two different constitut¬ 
ing magnetic ions. While the quasi-plateau phenomenon 
should still remain a rather subtle effect in heterodinu- 
clear complexes composed of Cu^+ (spin- 1 / 2 ) and Ni^+ 
(spin- 1 ) magnetic ions due to a relatively small difference 
between the g-factors not exceeding a few percent it 
should become much more pronounced in heterodinuclear 
complexes composed of Co^+ (spin- 1 / 2 ) and Ni^+ (spin- 
1 ) magnetic ions having much greater difference between 
g-factors (typically gco ~ 5.9 and gm ~ 2 . 3 )i ^^i^^i^^ 


V. SPIN-1/2 HEISENBERG TRIMER 


The next by complicity spin system with the different 
(/-factors is the triangle with uniform coupling and with 
only two (/-factors given by the Hamiltonian 


nt„m = JiS^S^ + SfS^ + ASfSI ( 31 ) 
-b S^S^ + SfSy + AS^SI 
+ + AS^,SI) 

- giBSt - 92B (S^ + SI). 


l'i^i) = int), 1 ^ 2 ) = I ;u), (34) 

|4'3) = |t)l|^)23, 14 / 4 ) = I ;)l 15)23, 


|4'7,8) = —j= 

\J‘^ + £± 

where the number in the lower right angle of the sym¬ 
bol I) correspond to the certain spin in the triangle, and 
\S), |T±) and jTo) are spin-singlet and components of the 
spin-triplet: 


(721 t)l|Ti))23+C±U)l|r+)23) , 
( 72 | ;)i|To)23 + c±|t)i|r_)23), 


I5) = ^(in)-ut)), 

\T+) = in), 

\T-) = I n), 

iTo) = ^ (I n) +1 n)). 

And the coefficients are 

-JA + 2 g.B±Q+ 


c± = 

c± = 


JA 

-JA- 2 g-B±Q. 

JA 


( 35 ) 


( 36 ) 


For our purposes the eigenvectors from j'l's) to I^Fg) are 
of special interest, as they demonstrate the monotonous 
explicit dependence of the magnetization on the magnetic 
filed under constant value of the projection of the total 
spin, which is in our case. The expectation value of 
the magnetic moment for the eigenstate with the lowest 
energy among the others (Id's) for positive B) is 


( 4 / 8 |{gi 5 J + (/ 2 ( 5 | + 5 |)}|d/ 8 ) 
1 2 gi - (51 - 252 ) 7 
“2 2 + ci 


( 37 ) 
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It is easy to see that for the case of the uniform g-factors, 
91 = 92 = 9 the magnetic moment expectation value be¬ 
come a constant equal to g/ 2 . Though, for the eigenval¬ 
ues of the magnetization operator the limit g- = 0 gives 
the correct result, this is not the case for the eigenvec¬ 
tors of the Hamiltonian. Thus, one cannot obtain the 
standard basis for the spin trimer by putting gi = g2 
in the Eqs, (IMl) . The susceptibility for the continuous 
magnetization dSzl) is given by 

+ 52 (^ 2 ^ + ^ 1 )} Id/g) = ( 38 ) 

45-c_ dc^ 

^ ~ {2 + cl)^~dB' 


The susceptibility becomes zero at gi =52- The partition 
function for the spin-trimer under consideration can be 
obtained in a straightforward way: 


= 2 ePi 






2 J 


e-^^ch{l3Q+) + e^^ch{l3Q-) 


( 39 ) 


The finite-temperature magnetization reads 


J 


Mt, 


0 - 
gp 4 


g p 4 


(ffi +252)6 ^-^sh +^92)B '\j +gie^'’^ sh 


( 40 ) 


5 S 2 B 


52 ch(/ 3 Q_)- 


g.{J- 2 g_B) 

Q- 


sh(/ 3 Q_) 


- e (52ch(^Q+)- ^ ^^ sh(/ 3 Q+) 


The plots of the normalized zero-temperature magneti¬ 
zation are presented in Fig. 0 Here the development 
of the non-plateau part of the magnetization curve with 
the increase in the difference 52 — 51 is clearly visible. 
For comparison the ordinary curve for 51 = 52 is also 
presented with a plateau at M/Msat = 1 / 3 , which corre¬ 
sponds to the ground state with = 3 / 4 , = 1 / 2 : 

^ (I m) +1 nt) +1 ;n)), ( 41 ) 

which transforms into the j^'g) 92 9 i- Let us men¬ 

tion also, that the B = 0 ground state at for 51 = 52 
is four-fold degenerate and the magnetic field lifts this 
degeneration just partly, because the 1/3 plateau state 
is still two-fold degenerate. The deviation from the hor¬ 
izontal line becomes more pronounced with the growing 
difference between g-factors of the spins. Thus, the zero- 
temperature magnetization curve for the simple system 
with finite discrete spectrum mimics the magnetic be¬ 
havior of magnets with the band of magnetic excitations. 
The value of critical field at which the level crossing be¬ 


tween iTg) and the fully polarized state [Ti) takes place 
is 


12 _ j‘^9+^ — 91 + 

Jdq — J 


\/(25_A -51)2 -H 85152 
45152 


( 42 ) 


In the limit 51 = 52 = 5 the value of critical field is 
The interplay between thermal fluctuations and 
the non-plateau behavior can be seen in the Fig. [ 5 ]where 
one can see a gradual smearing out of the magnetization 
curve with the rise of the temperature. 


VI. ISING-HEISENBERG (ISING-XYZ) 
DIAMOND CHAIN WITH DIFEERENT LANDE 
G-FACTORS 


To illustrate the features of having the spin cluster 
with non-conserving magnetization as a constituent of 
the Ising-Heisenberg one-dimensional systems let us now 
consider the simplest Ising-Heisenberg spin-chain with 
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- S2 = ^-^ - - - ^ 2 ^ 


1 


0.8 


M 0.6 
M 

sat 

0.4 


0.2 


0 

0 0.5 1 1.5 2 

B/J 

FIG. 7: The Zero temperature normalized magnetization 
curves for the Heisenberg spin trimer with two different Lande 
g-factors for J = 1, A = 2, = 2 and g 2 = 2 (solid red); 

g 2 = 4 (dotted black); 32 = 6 (dashed blue) and 32 = 10 
(dot-dashed brown). Msat = ^ (fli + 252 ). 



T 

T 

T 

T 


— = 0.5 - 

- — = 0.1 — 

— = 0.00001 

J 

J 

J 

J 



FIG. 8 : Finite temperature normalized magnetization curve 
for the Heisenberg spin triplet with two different Lande g- 
factors for J = 1, A = 2, gi = 2, 32 = 4 and T = 1.4J (solid 
red); T — 0.5J (dotted black); T = O.IJ (dashed blue) and 
T = O.OOOOIJ (dot-dashed brown). Msat = |((?i + 232 ) = 5. 


the XYZ-dimers and different g-factors, the diamond 
chainjii^— However, we are not going to describe the 
whole problem in all details, this can be a topic of the 
forthcoming and separate investigation. We just want 
to illustrate how rich the structure of the magnetization 
curve can be, if we include the spin cluster with non- 
conserved magnetization into the more involved struc¬ 
tures. The interest toward the diamond chain is large not 
only because of the simplicity of the system, especially 
in case of Ising-Heisenberg one-dimensional systems, but 



► ( 

f S„,i 

1 i 

• 





/ \ 

/ ^ n+\ 




> S„a > 



FIG. 9: The Ising-Heisenberg diamond-chain. Solid lines rep¬ 
resent the quantum interactions, while the dotted once stand 
for the interaction involving only ^-components of the spins. 
Here we also consider the y-factors of the quantum spins 
and Sj ^2 to be gi and (j 2 respectively. 


also as the diamond chain is believed to be the real mag¬ 
netic structure of the mineral azurite 42 r— The lattice 
in depicted in Fig. | 9 l where the quantum spin-dimer 
are the vertical bonds(solid lines) while the dashed lines 
correspond to Ising couplings. The Hamiltonian for the 
whole chain is the sum over the block Hamiltonians: 

'tidc=J 2 \ ~ J ^ 9 j+i^j+i) > , ( 43 ) 

i=i ^ 

n, =J {{1 + 7) + (1 - 7) SISI2 + ASI.SI^} 

+ - 9 iB)Sli 

+ + ^j+i) ~ 92B)Sj 2, 

where the g-factors of the Ising intermediate spins alter¬ 
nating with the spin-dimer are also taken alternating 


(93, j is odd 
[94, j is even. 


( 44 ) 


A. Exact solution 


One can apply the standard technique of the 
generalized classical transfer matrix to calculate 
the free energy of the model under consideration 
exactlyHowever, as here we deal with 
the alternation of two kind of blocks one has to compose 
a two block transfer matrix just by multiplying transfer 
matrices for odd and even blocks. In the other words, 
the partition function of the model can be factorized in 
the following form (the cyclic boundary conditions are 
assumed): 

N 

Zdc = ^ Xrj 

cri,...crN i=l 

= E ( 45 ) 


where with the aid of the block Hamiltonian eigenvalues 
from Eq. ®, one can obtain four quantities 

= ± 1 / 2 , as the entries of the following matrix: 


V = 2 e 


JA 

4 


_£C 93+£42 7-T- ^(93-34) 

i U- eP i Uo 

a -r T o ^(S3+S4) 

P 4 Uo e P 4 


, ( 46 ) 


e 
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where 

U± =W + ch (^^^{Bg+± 2 Kf + , 

C/O =W + ch + J2^2^ ^ 

W ch . ( 47 ) 


Thus, the partition function can be rewritten in the 
form 


^ — / 5 - 

^dr. — fc- 


Tr Tt 


( 48 ) 


where, 


J 


T = 


O 5(83 + 84) ,.o Q -5(33 ^84) ,-n 

2 t/i + 2 Uq 

+ e-^^C/+} t/o 


+e-^^[/+}c/o \ 


Then, for the free energy per unit cell in the thermo¬ 
dynamic limit, iV —>• oo, we have, 

/ = ^ - ^ (log 4 -b log A), ( 49 ) 

where the A is the largest eigenvalue of the matrix T, 
which is expressed by the entries of the matrix T in the 
following form: 



Now, we will analyze the non-plateau magnetization 
for Ising-XYZ diamond chain with different g-factors. 
Using the free energy, one gets the magnetization 


M = - 



^1 

^A \^) 


(51) 


As the lattice has six spins in the translational invariant 
unit cell, two a spins and two vertical dimers, the total 
saturation magnetization per one block (note that N is 
the number of block which is supposed to be even, while 
the number of the unit cell with six spins is N/ 2 ) is 

Msat = 2 (^9i + 52 + 2 (ff3 + 54)^ • (52) 


One can see the plots of the zero-temperature magneti¬ 
zation of the Ising-XYZ diamond chain in Fig. (I 10 |) . The 
zero-temperature curves can be obtained as a sufficiently 
low-temperature limit of the exact expression m- One 


can see series of quasi-plateau and magnetization jumps 
which can undergo a certain variation under the changing 
of the XY-anisotropy 7 (upper panel) or axial anisotropy 
A (lower panel). 

B. Ground states 


We are not going to present the comprehensive anal¬ 
ysis of all possible ground states and all types of mag¬ 
netization curves for the XYZ-Ising diamond-chain with 
different g-factors, but let us just illustrate the ground 
state structure of the T = 0 magnetization curve from 
Fig. (flOh l To be specific let us consider the magnetiza¬ 
tion curve for J = l,iF = l,gi = 8,52 = 2,53 = 2,54 = 
4,7 = 0.5 and A = 1 . 2 . The comparison of the different 
combinations of the spin-dimer eigenstates (Eq. ([ 9 ]) with 
Dz = 0 ) with the orientation od the intermediate Ising 
spins gives rise to series of possible eigenstates for the 
chain. However, only few of them are realized in the case 
we considered here. First of all, the B = 0 ground state 
is macroscopically degenerate for our choice of parame¬ 
ters. In this ground state all vertical quantum dimers 
are in |'I'2) state form Eq. (jH]) and all intermediate Ising 
spins can freely point either up or down. Thus, we have 
here 2 ^ configuration with the same energy. Arbitrary 
but non-zero magnetic field lifts this degeneracy. The 
system passes through the following ground state with 
the increasing the magnetic field which correspond to the 
course of the magnetization curve we are analyzing here: 


IGS”!) ^ |G 5 ' 2 ) ^ IGS”!) ^ \GS 3 ) \QS), ( 53 ) 


where the eigenstates with the corresponding magnetiza¬ 
tion and energies are 
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N 


i=i 

Bgl 


( 54 ) 


Ml = 


JA 1 


I - + ^ (53 + 54) , El = + J^-Yi93+ 54) , 

2 ^ B^gl + .P 4 4 2V 4 




N 


iG52)=n I i). ® 

i=i 


M2 = 


i^i^^ 2 ^^=-\{g 3 +g 4 ), E2 = ^- yiBg+ + 2 Kf + J272 + I (33 + 34 ), 


2 B^{Bg+ + 2 Kf + Pj 


N/2 

|G 53 ) = H I i)2j-l 0 | 4 ' 4 ) 2 j -1 0 I t)2j 0 | 4 ' 4 ) 2 „ 
i=i 


7 W 3 = 


S5 


JA 1 


2 ^/B^gl + J272 


- 7 (j3 - J4) , E3 = — - -JB^gl + J272 + _ (^3 - g^) 


B 


N 


IQ5)=n 11), 0 

i=i 

(Bg+ - 2 Kf 


JA 1 


■^QS = - , ' + 7 (53 + 54) , Gqs = — - {Bg+ - 2 Kf + y (53 + 34), 

2 B^{Bg+- 2 Kf + J^^^ ^ ^ ^ ^ 


here, 14^2,4)4 stand for the corresponding eigenvectors of taking into account the influence of the interaction with 
the spin-dimer (Eq. ([ 9 |) with = 0 ) at j-th block and the neighboring Ising spins aj and cr^+i, which leads only 
the | 4 'J)j is structurally the same |'t'4)j eigenstate but to a modification of the coefficient B_: 


^ ^ (I tt) + Bt\ U)), ( 55 ) 

^ Bg+ ± 2 K- ^{Bg+± 2 Kf + 

B_ = ---, 


the arrows indicate the orientation of the corresponding 
Ising spins. Here several comments are in order. First 
of all, one can see the phenomenon of reentrant phase 
transition when the system with increasing the magni¬ 
tude of the magnetic field enters the same ground state, 
IGS"!), twice. It also clearly can be seen from the ground 
state phase diagram which is presented in Fig. CB ■ 
The return to the ground state [G^l) is possible due to 
non-linear magnetic field dependence of the correspond¬ 
ing magnetization and the energies of all ground states. 
One can also see from the Eq. (l 54 )) that despite the visible 
ideal horizontal character of some part of the magnetiza¬ 
tion curve, they are not the magnetization plateaus, but 


just the very slowly growing part of the curve. Thus, the 
magnetization always has an explicit dependence on the 
magnetic field, unless 7 = 0 or/and gi =52- The quasi- 
saturated state, I< 3 /S'), at finite values of 7 has the mag¬ 
netization asymptotically converging to Mgat- However, 
this value is inaccessible for finite values of magnetic field. 
Although most of quasi-plateaus appear below the satu¬ 
ration (last critical) field, there also exists an alternative 
mechanism for quasi-plateau formation. Namely, the last 
plateau emerging above the last critical field (which is in 
fact not a true saturation field) may change to the quasi¬ 
plateau due to the XY anisotropy and consequently, the 
magnetization varies continuously above the last criti- 
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cal field and it never reaches full saturation except for 
asymptotically infinite magnetic field. 


-■-A = 0.9 - A =1 -A=l.l - -A= 1.2 



-y= 0. - y= 0.4 - y= 0-8 - Y= 10 



B/J 


FIG. 10: Zero temperature normalized magnetization as a 
function of magnetic field. Assuming fixed values: J = 1, gi = 
8 , p 2 = 2 ,53 = 2, g 4 = 4, AT = 1 (a) For a range values of A 
and fixed 7 = 0.5. (b) For a range values of 7 and fixed 

A = 1.2. Msat = f (51 + 52 + ^ (93 + 94)) = 13/2. 


VII. CONCLUSION 

In the present work, we have investigated in detail 
an absence of actual plateaus in zero-temperature mag¬ 
netization curves of quantum spin clusters and chains. 
It has been convincingly evidenced that strict plateaus 
may disappear from a magnetization process on assump¬ 
tion that constituent spins of quantum spin clusters have 
different Lande g-factors. More specifically, we have 



FIG. 11: Zero temperature groud-states phase diagram A 
versus B/J of the XYZ-Ising diamond-chain withe different 
g-factors for A = 1, 7 = 0.5 and g-factors 51 = 8, 92 = 2, 
93 = 2 and 54 = 4. 


demonstrated this intriguing feature on a few paradig¬ 
matic examples of quantum spin clusters such as the spin- 
1/2 Heisenberg dimer, the mixed spin-(l/ 2 ,l) Heisenberg 
dimer and the spin- 1/2 Heisenberg trimer, whereas the 
same phenomenon has been also found in the spin- 1/2 
Ising-Heisenberg diamond chain. From this perspective, 
the absence of actual magnetization plateaus due to the 
difference in Lande g-factors can be regarded as a gen¬ 
eral feature of low-dimensional quantum antiferromag- 
nets, because it emerges whenever the total magnetiza¬ 
tion does not represent a conserved quantity with well 
defined quantum spin numbers (the total magnetization 
need not be proportional to the total spin). Accord¬ 
ingly, a smooth variation of total magnetization can be 
attributed to a nonlinear dependence of a few (or all) 
discrete energy levels on a magnetic field. 

A few remarks are in order here, which might be useful 
for possible experimental testing of this interesting phe¬ 
nomenon. Although the magnetization curves of quan¬ 
tum spin clusters without strict plateaus may mimic to a 
certain extent the magnetization curves of quantum spin 
chains with continuous energy bands, it is obvious that 
the absence of magnetization plateaus does not in turn 
mean a gapless excitation spectrum. On the contrary, 
small magnetic spin clusters should always have an en¬ 
ergy gap, which can be easily experimentally tested by 
various resonance techniques. 

It is also worth noting that the deviation of magnetiza¬ 
tion from a strict plateau is proportional to a difference 
between the Lande g-factors, which makes an experimen¬ 
tal verification of this phenomenon more difficult. As a 
matter of fact, the most of transition-metal ions as for in¬ 
stance Cu^’*', Ni^+, Mn^+, Cr^+ or high-spin Fe^'*' with 
zero or totally quenched angular momentum can be de¬ 
scribed by the notion of more or less isotropic quantum 
Heisenberg spins and hence, these magnetic ions usu¬ 
ally have g-factors quite close to the free electron value 
g Tn 2 ^— Under these circumstances, it is custom¬ 
ary to combine the local Lande g-factors of individual 
magnetic ions into a single molecular g-factor as long as 
the isotropic Heisenberg exchange significantly prevails 
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over the zero-field splitting, asymmetric and/or antisym¬ 
metric exchangei^ Many experimental studies focused 
on a magnetism of such oligonuclear complexes there¬ 
fore simply ignore different Lande g-factors of individ¬ 
ual magnetic ions as the isotropic exchange is by far 
the most dominant coupling. However, a few transition- 
metal ions with unquenched angular momentum may 
have much higher Lande g-factors due to a relatively 
strong spin-orbit coupling like for example the low-spin 
Fe^’*' ion with typical value of g ~ 2.8 or Co^"*" ion with 
g Ki 6.0i^“— Another possibility how to increase the dif¬ 
ference of the Lande g-factors in oligonuclear complexes is 
to combine the almost isotropic transition-metal ion with 
highly anisotropic rare-earth ions, which may even pos¬ 
sess much greater Lande g-factors (e.g. Dy^+ typically 
has g Ri 20) though this extraordinary large g-value usu¬ 
ally correlates with a rather strong anisotropy in the ex¬ 
change interactioui ^^i^^i^^ Under the extreme situation, 
the XY-part of exchange coupling might be even of op¬ 
posite sign than the Z-part (ferromagnetic versus anti¬ 
ferromagnetic) as it has been recently reported for the 
heterodinuclear Cr^+-Yb^+ complexi^ 

Last but not least, let us briefly comment on experi¬ 
mental implications for the quantum spin clusters studied 
in the present work. The spin-1/2 Heisenberg dimer has 
previously proved its usefulness as the plausible model 
of many homodinuclear Cu^+-Cu^+ complexes 
However, the difference between the local Lande g-factors 
in the homodinuclear coordination compounds may only 
stem from a different coordination sphere of individ¬ 
ual magnetic ions, whereas this difference does not ex¬ 
ceed in most cases a few percent that would be insuf¬ 
ficient for an experimental testing. On the contrary, 
the considerable difference in the local g-factors could 
be expected in heterobimetallic coordination compounds, 
which are composed of the nearly isotropic magnetic ion 
(e.g. Cu^+, Ni^+, high-spin Fe^+, etc.) and the highly 
anisotropic magnetic ion (e.g. Co^’*', low-spin Fe^+, etc.). 
Hence, the heterodinuclear Co^+-Cu^+ and Fe^“''-Cu^+ 


complexes could be regarded as sought experimental re¬ 
alizations of the generalized spin-1/2 Heisenberg dimer, 
which may display a substantial deviation of the mag¬ 
netization from a strict plateau as the g-factors of in¬ 
dividual magnetic ions may even possess opposite signs 
due to a spin-orbit coupling (e.g. gpe ~ ~l-7, gcu ~ 2.1 
was reported in Ref. , and negative g-factors of Co^'*' 
and Cu^’*' magnetic ions were investigated in Ref. [13). 
The similar conjecture can be formulated for experimen¬ 
tal representatives of the mixed spin-(1/2,1) Heisenberg 
dimer. As usual, the magnetic anisotropy in heterodin¬ 
uclear Cu^“''-Ni^+ complexes does not cause a significant 
deviation of the magnetization from a strict plateauj^i^ 
but a rather large deviation could be expected instead 
in heterodinuclear Co^'''-Ni^+ coordination compounds 
with typical values of the g-factors gco ~ 5.9 and gNi ~ 
2.3.— It can be also anticipated that the homotrinuclear 
CiU+-Cu^+-CiU+ complexes^^— as experimental repre¬ 
sentatives of the spin-1/2 Heisenberg trimer should not 
have a significant deviation of the magnetization from 
a strict plateau unlike the heterotrinuclear Cu^^-Co^"*"- 
Cu^'*' complexes— 
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